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(57) Abstract 



The invention relates to a method and system for 
carrying out stock management in a multilevel distri- 
bution chain by applying optimisation algorithms to a 
range profile built on a distribution network formulation 
and on different elements which are invoiced in the net- 
work. The optimisation problem is defined in a formal 
way and algorithms for maximising stock autonomy by 
means of requirement profiles are proposed. Said op- 
timisation process allows transport costs to be reduced 
to a minimum and bottlenecks to be avoided in stocks. 
But should bottleneck occur, the system allocates prior- 
ities according to need and provides an optimum solu- 
tion for monotonous product distribution. The inventive 
method comprises the following steps: initialisation of 
stock autonomy profile with a starting value; calcula- 
tion of required entries to satisfy profile restrictions; 
creation of the most cost-effective flow for said needs 
by means of a minimal autonomy of the requirement 
profile until a solution is found. 
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(57) 7jreanimen fagging 

Verfahren imd System z.um Vcrwaltcn von Bestftnden in ciner Mehifachebenendistributionskctte durch Anwcnden von Optimierungsal- 
gorithmen auf ein aus einer Foimulieiung dcs Distributionsnetzwerks und den in das Netzwerk faktoriertcn verschiedenen Elementcn kon- 
stmiertes Bereichsprofils. Eine formale Definition des Optimicrungsproblcms wird definiert, und Algorithmen fQr die Maximiemng der 
Lagerreichweite yon Abdeckungsprofilen werden vorgeschlagen. Durch den OptimienmgsprozeB werden Traxisportfcosten auf ein Minimum 
reduzieit und Bestandsengpasse vennieden. Falls jedoch Engpflsse auftreten sollten, veisieht das System Bedarfe mit Priorit&en und liefert 
eine opu'male Losung fQr eine monotone Verteihing von Produkten. Das Verfahren weist folgende Schritte auf: Initialisieien des Reich- 
weitenprottls des Lagers mit einem Startwert; Berechnen der notwendigen Zuflusse, die notwendig sind, urn die ProtUbeschrankungen zu 
erfullen; Konstruieren des kostengOnstigsten Fhisses fur diese Bedarfe unter Verwendung eines Minimalkostenalgorilhmus; und Verieleinems 
oder VergroBems der Lagerreichweite der Abdeckungsproflle, bis eine Losung gefunden ist. 



LEDIGUCH ZUR INFORMATION 

Codes zur Identifizierung von PCT-Vertragsstaaten auf den Kopfbogen der Schriften, die internationale Anmeldungen gemass dem 
PCT veroffentHchen. 



AL 


Albanien 


ES 


Span ten 


LS 


Lesotho 


SI 


Slowenien 


AM 


Annenien 


Fl 


Fmnland 


LT 


Litauen 


SK 


Slowakei 


AT 


Osterreich 


FR 


Frankreich 


LU 


Luxemburg 


SN 


Senegal 


AU 


Australien 


GA 


Gabun 


LV 


Lett land 


SZ 


Swasiland 


AZ 


Asetbaidschan 


GB 


Vercintgtes Kflnigreich 


MC 


Monaco 


TD 


Tschad 


BA 


Bosnien-Herzegowina 


GE 


Georgien 


MD 


Republik Moldau 


TG 


Togo 


BB 


Barbados 


GH 


Ghana 


MG 


Madagaskar 


TJ 


TadschikUtan 


BE 


Belgien 


GN 


Guinea 


MK 


Die ehemaltge jugoalawische 


TM 


Turkmenistan 


BF 


Burkina Faso 


GR 


Griechenland 




Republik Mazedonien 


TR 


TOrkei 


BG 


Bulgarien 


HU 


Ungam 


ML 


Mali 


TT 


Trinidad und Tobago 


BJ 


Benin 


IE 


Irland 


MN 


Mongolei 


UA 


Ukraine 


BR 


Brasilien 


IL 


Israel 


MR 


Mauretanien 


UG 


Uganda 


BY 


Belarus 


IS 


bland 


MW 


Malawi 


US 


Vereinigte Staaten von 


CA 


Kanada 


IT 


II alien 


MX 


Mexiko 




Amerika 


CF 


Zeruralafrikanische Republik 


JP 


Japan 


NE 


Niger 


uz 


Usbekistan 


CG 


Kongo 


KE 


Kenia 


NL 


Niederlande 


VN 


Vietnam 


CH 


Schweiz 


KG 


Kirgisistan 


NO 


Norwegen 


YU 


Jagoslawien 


CI 


Cote d'lvoire 


KP 


Demokraiische Volks republik 


NZ 


Neuseeland 


ZW 


Zimbabwe 


CM 


Karri erun 




Korea 


PL 


Polen 






CN 


China 


KR 


Republik Korea 


PT 


Portugal 






CO 


Kuba 


KZ 


Kasachstan 


RO 


RumSnien 






cz 


Tschechische Republik 


LC 


Sl Luda 


RU 


Russische Federation 






DE 


Deutschland 


U 


Liechtenstein 


SD 


Sudan 






DK 


DSnemark 


LK 


Sri Lanka 


SE 


Schweden 






EE 


Estland 


LR 


Liberia 


SG 


Singapur 








WO 99/67729 PCT/EP99/04229 



Verfahren und System zur Maximi erung des Bereichs von Ab- 
deckungsprof ilen bei der Bestandsverwaltung 

S 

Die vorliegende Erfindung betrifft ein Verfahren zur Opti- 
mierung der Transportplanung und insbesondere ein Verfah- 
ren zur Optimierung der Lagerreichweiten in einem mehrere 
10 Verteilebenen aufweisenden Transportnetzwerk durch die 
Verwendung eines rechnergestutzten Optimierungsalgorith- 
mus, insbesondere zur schnellen Bestimmung von Terminpla- 
nen fiir Versand und Produktion. 

15 Allgemeiner Stand der Technik: 

In vielen Indus triebereichen, die die Fertigung und Ver- 
teilung von Gutern betreffen, besteht ein Bedarf an der 
Zuteilung von Bestands- und Transportressourcen . Diese Zu- 

20 teilungsentscheidungen sind in der . Regel Beschrankungen 
hinsichtlich Ausrustung, Zeit, Kosten, Lagerkapazitat und 
anderen Parametern unterworfen, die das Ergebnis eines Di- 
stributionsprozesses, insbesondere eines Warenver tei lungs - 
prozesses, beeinf lussen. Als Beispiel, das hier von beson- 

25 derem Interesse ist, besteht ein Bedarf daran, die Distri- 
butions- und BestandsgroEen einer Versorgungskette zu op- 
timieren, in der mehrere Produktionsstatten existieren, 
wobei mehrere Verteilzentren strategisch in der Nahe von 
Kunden lokalisiert sind. 

30 Entscheidungen uber die Zuteilungen von Ressourcen 

sind in der Regel Beschrankungen unterworfen. Einerseits 
sind Ressourcen von ihrer Gesamtverfugbarkeit immer be- 
grenzt, und aufierdem kann die Niitzlichkeit einer bestimm- 
ten Ressource in einer bestimmten Anwendung ebenfalls be- 

35 grenzt sein. So ist beispielsweise die Verkehrsfuhrungs- 
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kapazitat jeder einzelnen Verbindungsstrecke in einem Te- 
lekommunikationssystem begrenzt, w&hrend der dem Kommuni- 
kationssystem angebotene Gesamtverkehr ebenfalls begrenzt 
ist. Die einzelne Zuteilung einer Ressource weist eine 
5 "Rendite" auf, d.h. die Kosten oder den Gewinn fur diese 
Zuteilung. Das Problem besteht somit darin, die Ressourcen 
derart zuzuteilen, dafi alle Beschrankungen berucksichtigt 
sind und gleichzeitig die Rendite maximiert wird, d.h. die 
Nachteile, z. B. Kosten, minimiert bzw. die Vorteile, z. 

10 B. Gewinne, maximiert werden. 

Ein Verfahren, derartige Zuteilungsentscheidungspro- 
bleme darzustellen, ist als das Lineare Programmiermodell 
bekannt. Ein derartiges Modell besteht aus einer Reihe von 
linearen Beziehungen, die in einer Matrixform dargestellt 

15 werden und quantitativ die Beziehungen zwischen Zuteilun- 
gen, Beschrankungen und den Ergebnissen des Optimierungs- 
prozesses reprasentieren. In den linearen Beziehungen wird 
die Summe aus konstanten Koef f izienten multipliziert mit 
unbekannten Zuteilungswerten dargestellt. Eine derartige 

20 Modellierung der Linearen Programmierung wird in einem 
mehrdimensionalen Raum erzielt, wobei mehrdimensionale 
Vektoren eine mehrdimensionale Figur bzw. ein Polytop lie- 
fern, wobei jede Facet te auf einer Flache davon durch 
Gleichungen begrenzt wird, die Beziehungen zwischen zuge- 

25 teilten Ressourcen in dem Prozefc definieren. Ein Losungs- 
ansatz des linearen Programmierproblems stellt beispiels- 
weise der Simplexalgorithmus dar, der 1947 von George 
Dantzig entwickelt wurde. Eine Losungsalternative ist der 
Karmarkar-Algorithmus . Es gibt weitere Optimierungsalgo- 

30 rithmen, mit denen Mindestkosten-Flufinetzwerkprobleme ge- 
lost werden konnen. Netzwerkf lufipr obi erne stellen dabei ei- 
nen Sonderfall der linearen Programmierung dar. 

Netzwerkf lufcalgorithmen finden bei Planungsproblemen 
in der Industrie vielseitig Anwendung. Sie k6nnen bei- 
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spielsweise bei Problemstellungen hinsichtlich Zuordnung, 
Transport, Mindestkostenf lufc, kurzestem Weg und grofitem 
Flufc durch ein Netzwerk angewandt werden. 

Bei Zuordnungsproblemen gibt es einen zweiseitigen 
(bipartiten) Graph, der eine Anzahl produktiver Knoten und 
verbrauchender Knoten (beispielsweise Arbeitern und Aufga- 
ben) umfa£t. Die Bogen zwischen den Arbeitern und den Auf- 
gaben werden mit Kosten -fur die Zuordnung des Arbeiters zu 
der Aufgabe gewichtet. Das Optimierungsproblem besteht bei 
diesem Szenarium in der Zuordnung einer Aufgabe zu jedem 
Arbeiter derart, daS die Gesamtkosten auf ein Minimum re- 
duziert werden. Bei Transportproblemen gibt es einen zwei- 
seitigen Graph, bei dem die produktiven und verbrauchenden 
Knoten Fabrikanlagen und Verteilzentren sind. Jede Fabrik- 
anlage produziert Einheiten, und jedes Distributionszen- 
trum hat einen Bedarf an diesen Einheiten. Die Bogen zwi- 
schen Fabrikanlagen und den Verteilzentren werden durch 
Transportkosten der Einheiten gewichtet. Das Optimierungs- 
problem besteht in diesem Fall darin, die Transportkosten 
auf ein Minimum zu reduzieren, unter der Auflage dafc alle 
Bedarf e in den Verteilzentren erf till t werden. Das Mindest- 
kostenf lufcproblem ist lediglich das Transportproblem mit 
Zwischenknoten in dem Netzwerk. Aufcerdem konnen die Bogen 
eine kleinste und grdfcte Kapazitat aufweisen. 

Bei dem Problem des kurzesten Weges liegt ein Graph 
mit positiv gewichteten Bogen vor. Das Optimierungsproblem 
besteht darin, zwischen zwei gegebenen Knoten in dem 
Transportnetzwerk den kurzesten Weg zu finden. Das Problem 
des grofcten Flusses durch ein Netzwerk ist Shnlich dem 
Transportproblem, mit der Ausnahme, da£ die Bogen zwischen 
den Knoten begrenzte Transportkapazitaten, aber keine 
Transportkosten aufweisen. Das Optimierungsproblem besteht 
darin, ohne Verletzung der Transportkapazitaten einen 
grofcten Flufi durch das Netzwerk transportiert zu bekommen. 
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Jedoch sind im Stand der Technik Netzwerkf lufcprobleme 
auf Netzwerke mit Knoten und B6gen dazwischen begrenzt, 
wobei der FluS durch Transportkosten zwischen ihnen mit 
einer linearen Strafe belegt wird. Es besteht ein Bedarf 
5 an einer Netzwerkf lufcldsung, bei der die Minimierung der 
Transportkosten eine niedrigere Prioritat hat und bei der 
die Maximierung des Bereichs von Abdeckungsprof ilen in den 
Verbrauchsknoten von herausragender Wichtigkeit ist. Die- 
ses Problem kann weder, wie im Stand der Technik, als li- 
10 neares Optimierungsproblem noch als Netzwerkf lu&problem 
formuliert werden. 



Kurze Darstellung der Erfindung: 



Die Aufgabe der vorliegenden Erfindung besteht in der 
Schaffung eines neuen Algorithmus fur die Maximierung des 
Bereichs von Abdeckungsprof ilen bei dem Verteilungspro- 
blem, das sich bei Planungssystemen in der industriellen 

20 Produktion ergibt. Es werden mehrere Algorithmen vorge- 
schlagen, die auf dieses Problem anwendbar sein sollen. 
Die Entscheidung hinsichtlich der Frage, welcher Algorith- 
mus fur ein gegebenes Verteilungsproblem angemessen ist, 
hangt von der Grofce des Problems und der maximal akzepta- 

25 blen CPU-Zeit fur die Berechnung ab. Die vorliegende Er- 
findung schlagt eine neue Formulierung des Netzwerkf lufi- 
problems vor, die unterschiedliche Transportarten, Zeit- 
und Terminbeschrankungen, Bedarf sprioritfiten und festge- 
legte Produktionsf lusse berucksichtigt . Die algorithmische 

30 Funktion wird dann auf diese Formulierung angewendet. 

Die Aufgabe des Algorithmus besteht darin, die freien 
Variablen der Bereichsprof ilf ormulierung derart zu wahlen, 
da£ zuerst der Bereich von Abdeckungsprof ilen maximiert 
wird und danach die Transportkosten minimiert werden. Der 
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vorgeschlagene Algorithmic kann als grundlegenden Baustein 
jeden beliebigen Mindestkostenf luSalgorithmus verwenden. 

Ausfuhrliche Beschreibung der Erfindung: 



Die vorliegende Erfindung betrifft ein Planungssystem zum 
Optimieren des Flusses von Waren und Gutern in komplizier- 
ten Distributionsnetzwerken unter Verwendung eines neuen, 
effizienten Algorithmus . Zunachst wird die Formulierung 
10 des Planungssystems erortert, danach werden die auf diese 
Formulierung anzuwendenden neuen Algorithmen eingefiihrt. 
Beim Aufbau des Szenariums, in dem der Algorithmus ange- 
wendet werden soli, wird von der Annahme ausgegangen, daS 
das fur die betrachtete Zeitperiode zu verteilende Produk- 
15 tionsvolumen bereits festgelegt ist und da£ auch Bedarfs- 
daten bekannt sind. Der hier beschriebene Verteilungsalgo- 
rithmus soil hauptsachlich auf Storungen im Produktions- 
prozeE oder auf kurzfristige Veranderungen im Verteilungs- 
prozefi reagieren. Das System eignet sich gut fur Produkte, 
20 die in groBen Mengen produziert werden. und bei denen der 
Bedarf stark schwankt . Der hier dargelegte Optimierungs- 
prozeS versucht, durch Zukunf tsplanung (Push-Distribution) 
Transportkosten auf ein Minimum zu reduzieren und Engpasse 
zu vermeiden. Wenn es dennoch zu Engpassen kommt, wird an- 
25 gestrebt, optimal zu reagieren (Fair-Share-Distribution) , 
indem fur den Bedarf Prioritaten gesetzt werden. 

Die grundlegende Struktur eines Systems zum Verteilen 
von Waren, auf die die vorliegende Erfindung angewendet 
werden kann, umfafct aus drei Elemente: die Produktionsan- 
30 lagen, die zentralen Lagerhauser und die Verteilzentren. 
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Kurze Beschreibung der Zeichnungen: 

! ,:;;vDie Erfindung wird nachfolgend anhand in. Figuren darge- 
.. stellter Ausfuhrungsbeispiele naher erlautert. Es zeigen: 

5" 

Figur 1 eine graphische Darstellung eines einfachen Di- 
stributionsnetzwerkes . 

Figur 2 eine grapische Darstellung eines komplizierten Di- 
stributionsnetzwerkes . 
10 Figur 3 eine graphische Darstellung des Zuordnungsproblems 
in Form eines zweiseitigen Graphen. 
- Figur 4 eine graphische Darstellung eines Netzwerkes mit 
ausgeglichener Distribution. 

Figur 5 eine graphische Darstellung von Produktions- und 
15 Bedarf sknoten mit verschiedenen moglichen Zuordnungen. 
Figur 6 eine Produktionsanlage mit Angebot und Bedarf. 
Figur 7 eine graphische Darstellung eines Distributions- 
netzwerkes mit mehreren Transportarten . 

Figur 8 eine graphische Darstellung eines optimierten Di- 
20 stributionsnetzwerkes mit mehreren Transportarten. 

Figur 9 ein Flufcdiagramm sowohl moglicher als auch opti- 
mierter Transportstrecken. 

Figur 10 ein weiteres Beispiel moglicher und optimierter 
Transportstrecken . 

25 Figur 11 den Einflufc unzuiassiger Perioden auf die mogli- 
chen Transportoptionen in einem Distributionsnetzwerk. 
Figur 12 eine graphische Darstellung eines einfachen Di- 
stributionsnetzwerkes mit Zeitwerten auf den Bogen. 
Figur 13 eine graphische Darstellung der Produktions- und 

30 Bedarf sprofile eines Distributionsnetzwerkes . 

Figur 14 eine graphische Darstellung eines globalen Be- 
reichsprof ils fur eine mogliche Zuordnung. 
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Figur 15 eine graphische Darstellung fur ein globales Be- 

reichsprofil fur eine weitere mogliche Zuordnung. 

Figur 16 ein Beispiel fiir einen monotonen optimalen Trans- 

portplan. 

Figur 17 die Suche nach dem maximalen monotonen Bereich- 
sprof il . 

FIG. 1 zeigt den einfachsten Fall eines Distributi- 
onsnetzwerkes . In einer Produktionsanlage 102 konnen unter 
Einsatz von Produktionseinrichtungen, welche begrenzt ver- 
fiigbar sind Waren produziert werden. Diese Waren werden 
anschliefcend (z. B. in Lagerhausern) zwischengelagert 
104. Urn kurze Warentransportstrecken zu Kunden und um ei- 
ne schnelle und wirtschaf tliche Anlieferung zu garantie- 
ren, sind Zwischenlagereinrichtungen bzw. Verteilzentren 
106 erf orderlich. Mit den Zwischenlagereinrichtungen (Zwi- 
schenknoten) bzw. den zentralen Lagerhausern 104 kann ein 
groSeres Transportvolumen erzielt werden. Der Einsatz zen- 
traler LagerhSuser 104 fiihrt zu langeren Anlief erzeiten, 
aber auch zu niedrigeren Transportkosten. Ein Verteilzen- 
trum kann, falls erf orderlich, auch eine direkte Anliefe- 
rung erhalten. Diese besondere Darstellung gibt nur eine 
schematische Struktur wieder. In der Praxis konnen Netz- 
werkstrukturen viel grofiere Abmessungen mit einer groSen 
Anzahl von Anlagen, Lagerhausern und Verteilzentren errei- 
chen. Wie in FIG 2 dargestellt, umfassen diese Netzwerk- 
strukturen mehrere Herstellungsanlagen 202, zentrale La- 
gerhauser 204 und Verteilzentren 206. 

Die Aufgabe des Verteilungsproblems kann wie folgt 
zusammengefafit werden. In der typischen Versorgungskette 
gibt es Quellen, Senken und ein Transportnetzwerk. Die 
Quellen sind Anlagen, fur die ein geplantes Produktions- 
profil vorliegt. Die Senken sind die Verteilzentren, fur 
die ein Bedarf sprof il (Kundenauf trage oder ein prognosti- 
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zierter Bedarf) vorliegt. Das Transportnetzwerk besteht 
aus Bogen mit Transportkosten und Zeitwerten zwischen den 
Zwischenknoten. Die vorliegende Erfindung sucht nach der 
optimalen Zuordnung der raumlich und zeitlich verteilten 
5 Produktion zu dem ebenfalls raumlich und zeitlich verteil- 
ten Bedarf, so dafi die Lagerreichweite von Abdeckungspro- 
filen maximiert wird. Wie oben erwahnt, werden zwei Falle 
unterschieden. In dem "Push" -Fall ubersteigt die Anzahl 
der zur Verfugung stehenden Produkte den Bedarf, und die 

10 Produkte werden derartig verteilt, dafi der Bereich der Ab- 
deckung von jedem gleichformig zunimmt . In dem "Fair Sha- 
re" -Fall liegt die Anzahl der zur Verfugung stehenden Pro- 
dukte unter dem Bedarf, und die Produkte werden derartig 
verteilt, da£ wichtigere Kunden bevorzugt belief ert werden 

15 und Kunden mit gleicher Prioritat ahnlich behandelt wer- 
den. Die Erfindung verhindert, einen Kundenbedarf mit 
niedrigerer Prioritat zu bedienen, bevor ein Kundenbedarf 
mit hoherer Prioritat bedient wird. In dem "Push" -Fall 
wird das Minimum des Bereichs von Abdeckungsprof ilen auf 

20 alle Verteilzentren angewendet. Fur den "Fair Share"-Fall 
wird jede Priori tatsklasse sequentiell behandelt, und die 
Abdeckung wird maximiert, urn diese Bedarf e mit hoher Prio- 
ritat zu bedienen. Wenn daruber hinaus verschiedene Trans - 
portarten, Transportwege und Angebotsquellen zur Verfugung 

25 stehen, sollte weiterhin der gesamte Verteilplan so opti- 
miert werden, dafi die Gesamtkosten so niedrig wie moglich 
sind. 

Bei dem Ansatz zur L6sung des oben beschriebenen Ver- 
teilungsproblems spielt das Transportproblem bzw. das Zu- 
30 ordnungsproblem eine zentrale Rolle. Das zugrunde liegende 
Zuordnungsproblem ist wie folgt definiert: in einer zwei- 
seitigen graphischen Darstellung mit einer gleich groSen 
Anzahl linker und rechter Knoten sollte jeder Knoten auf 
der linken Seite einem Knoten auf der rechten Seite zuge- 
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ordnet sein, so da£ jeder Knoten genau einen Partner hat. 
Die resultierende Verbindung PBogen") zwischen zwei Kno- 
ten ist mit bestimmten Kosten verbunden. Die resultieren- 
den Gesamtkosten sollten minimal sein. FIG. 3 beschreibt 
das Zuordnungsproblem in Form einer graphischen Darstel- 
lung. Der Pfeil 302 stellt eine mogliche Zuordnung einer 
Quelle 304 auf der linken Seite zu einer Senke 306 auf der 
rechten Seite dar. Die Kosten von Zuordnungen sind als 
Werte 308 der B6gen eingetragen. Das klassische Transport- 
problem unterscheidet sich von dem Zuordnungsproblem da- 
hingehend, dafc ein Knoten auf der linken Seite (nun als 
"Produzent" bezeichnet) einer beliebigen Anzahl von Knoten 
auf der rechten Seite {nun als "Verbraucher" bezeichnet) 
zugeordnet sein kann. Das zu losende Problem besteht dar- 
in, wie zwischen den Knoten eine geeignete Verteilung von 
Materialien gefunden werden kann, so dafc der sogenannte 
Bestandsbereich maximiert wird. Dies bedeutet, den klein- 
sten Bestand uber alle Knoten und Perioden hinweg zu maxi- 
mieren. AuBerdem soil aus den moglichen Losungen (Losungs- 
schar) die Losung mit den geringsten Transportkosten ge- 
wahlt werden. Der Lagerbereich eines Lagerhauses ist als 
die Anzahl nacheinanderf olgender Perioden (einschliefclich 
der gegenwartigen Periode) definiert, fur die der aktuelle 
Bestand ausreicht. Ein Bereich von 1 bedeutet somit, dafc 
der Bestand exakt den Bedarf einer einzigen, namlich der 
gegenwartigen Periode abdeckt. 

Die Entwicklung des Problems, auf die die vorgeschla- 
genen Algorithmen angewendet werden sollten, wird nun na- 
her beschrieben. 



DER KLASSISCHE FALL 

Im einfachsten Fall betrachte man ein einfaches zweiseiti- 
ges Netzwerk mit Produktionsanlagen (PA) auf der einen 
Seite und Verteilzentren (Verteilzentrum) auf der anderen 
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Seite. In diesem klassichen Fall erfolgen die Anlieferun- 
gen immer von den Anlagen zu den Verteilzentren und nicht 
in der entgegengesetzten Richtung. Eine PA kann ein belie- 
biges Verteilzentrum versorgen. FIG. 4 zeigt ein Beispiel 
5 mit zwei Anlagen 402, 404 und zwei Zentren 406, 408. Nor- 
malerweise wird man wesentlich mehr Verteilzentren als 
Produktionsanlagen haben. 



t eine Produktionsrate pi(t) auf . In einem Verteilzentrum 
10 j liegt eine Kundennachfrage von Oj(t) vor. Der anfangli- 
che Bestand aller Anlagen und Zentren betragt Null. Wenn 
der Bereich fur 1 Periode gesetzt ist (d.h. der Bestand 
eines Verteilzentrums mu£ in einer bestimmten Periode nur 
genau fur diese Periode reichen) , kommt es zu folgender 
15 Situation: Der Bedarf eines Verteilzentrums j in einer Pe- 
riode t kann durch alle Produktionen abgedeckt werden, die 
rechtzeitig dem Verteilzentrum angeliefert werden konnen 
(das heiSt nicht spater als Periode t) . Urn diesen Bedarf 
abzudecken, konnte deshalb die Produktion aller Anlagen i 
20 in alien Perioden t ' < t-T±j (wobei T ±j die Transportzeit 
ist) verwendet werden. Bei diesem ersten Beispiel sind die 
Transportzeiten gleichformig auf eine Periode eingestellt. 
Beim Transport einer Materialeinheit entstehen Transport- 
kosten Cij . FIG. 5 zeigt in unserem Netzwerk diejenige Si- 
25 ' tuation, bei der das Fenster auf drei nachfolgende Peri- 
oden erweitert wird (wodurch in die Netzwerkberechnung die 
zusatzliche Variable der Zeit eingefuhrt wird) . Die Anzahl 
der beriicksichtigten Perioden wird als der Horizont be- 
zeichnet. Bei diesem neuen Szenarium steht jeder Knoten 
30 502 in der Darstellung nun fur eine Tagesproduktion bzw. 
einen Tagesbedarf. Die Pfeile 504 entsprechen nicht mehr 
direkt den Transportwegen, sondern entsprechen stattdessen 
der logischen Zuordnung zwischen einer bestimmten Produk- 
tionseinheit und einer bestimmten Bedarf seinheit . Jeder 



Eine Produktionsanlage i weist in einer Periode 
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Bogen kann deshalb durch den Ort der Quelle und des Ziels 
(i and j) und durch den Anfangs- und Endzeitpunkt (t* and 
tj) eindeutig identif iziert werden. Die Zuordnung einer 
bestimmten Materialmenge von einem Knot en auf der linken 
5 Seite zu einem Knoten auf der rechten Seite wird als der 
FlulS auf dem Bogen bezeichnet. Dieser Flufc fiihrt zu Kosten 
in der Hohe der entsprechenden Transportkosten C*j in dem 
Dis tributionsnetzwerk . 

Wie man sehen kann, ist dies ein klassisches, verein- 

10 fachtes Transportproblem, wie es aus Operations Research 
bekannt ist. Als solches kann es iiber einen standardmafii- 
gen Ansatz aus Operations Research gelost werden (bei- 
spielsweise liber den MODI-Algorithmus, wie er aus dem 
Stand der Technik bekannt ist) . Da der eigentliehe Opti- 

15 mierer auf der Ebene des Transportproblems arbeitet, wird 
das Transportnetzwerk von jetzt ab als das Optimierungs- 
netzwerk bezeichnet, und der Optimierer wird als der 
Transportoptimierer bezeichnet. Die Knoten in dem Verteil- 
netz werden als die Distributionsnetzwerkknoten bezeich- 

20 net, die Knoten in dem Optimierungsnetzwerk als Optimie- 
rungsnetzwerkknoten, wenn auf beide Bezug genommen wird. 
Die Knoten in dem Optimierungsnetzwerk werden in Quellen 
und Senken unterteilt. 

Dieser besondere Fall eines Distributionsnetzwerks 

25 ist einfacher zu behandeln, da hier Produktion und Bedarf 
des optimierenden Netzes genau denen in dem Distributions- 
netzwerk fur eine einzelne Periode entsprechen. Wenn der 
Bereich auf mehrere Zeitperioden vergrdfiert wird, ist dies 
nicht mehr der Fall. Die genaue Ubertragung dieser Werte 

30 aus dem Verteilnetz in das optimierende Netz wird in dem 
nachsten Abschnitt beschrieben. 
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VON DEM VERTEILUNGSPROBLEM ZU DEM TRANS PORTPROBLEM IN DEM 
' KLASSISCHEN FALL. 

In dem vorausgegangenen Abschnitt wurde das Vertei- 
lungsproblem fur einen festgelegten Bereich von einer Pe- 
riode auf das Transportproblem projiziert. Das Ziel be- 
steht nunmehr darin, die Projektion fur jeden moglichen 
Bereich zu verallgemeinern. Da die Werte der Senken in dem 
optimierenden Netzwerk bereits um einen Bereich von mehr 
als einer Periode von dem Bedarf in dem Verteilnetz abwei- 
chen und da, die Quellenwerte sich von den Produktionswer- 
ten in einem Distributionsnetzwerk unterscheiden konnen, 
wird eine weitere Unterscheidung hinsichtlich der Daten 
des optimierenden Netzwerks spezifisch von den Originalda- 
ten in dem Verteilnetz vorgenommen. Im Weiteren wird bezo- 
gen auf das optimierende Netzwerk immer von Angebot und 
Bedarf gesprochen - die Ausdriicke Produktion und Auftrag 
beziehen sich immer auf die Werte in dem Verteilnetz. Das 
Angebot in einer Periode t eines Knotens i in dem optimie- 
renden Netz ist als Si ( t) definiert, und der Bedarf ist 
als di(t) definiert. 

Definition des Bereichs der Abdeckung 

Der Zielbestand Li (t) in einer Periode t fiir einen gegebe- 
nen Bereich R(t), ergibt sich aus den Auftragen der nach- 
folgenden Perioden (die Einbeziehung des Sicherheitswaren- 
bestand mi (t) erfolgt hierbei noch nicht) : 

t+ lR(t)J 

Li(t)= I Oi(r) + (R(t) - LR(t)Joi(t-hfR(t) 7) 
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Alternative Definition / Verallgemeinerung des Bereichs 
der Abdeckung 

Die allgemeine Aussage uber den Bereich der Abdeckung (La- 
gerreichweite) kann alternativ so geschrieben werden, daS 
5 sie eine Maximierung des Multiplikators des Sicherheitswa- 
renbestands enthalt: 

Li(t) =Oi(t) +R(t) * Sicherheitswarenbestand 
oder allgemein fur jede monoton ansteigende Funktion f mit 

Li(t) = f (R(t) ) . 

10 

Der Ubersichtlichkeit halber wird im folgenden nur die 
erstgenannte Definition des Bereichs der Abdeckung be- 
trachtet . 

15 Im Fall eines Verteilzentrums hangt der Bedarf di(t) von 
der Differenz des Zielbestands Li(t) - Li(t-l) ab, der wie- 
derum von dem Bereich R(t)) und dem Kundenauf trag Ot(t)) 
abhangt. Der Bedarf di(t) ergibt sich zu: 

20 Verteilzentren: di(t) = Oi(t) + Lift;. - Li(t-l)' 

[Die Bedingungen: di(t)20 <=> R(t) > R(t-l) -1 sollten 
hierbei erfullt sein.) 

25 Das Angebot Si(t) einer Anlage entspricht bei diesem ein- 
fachen Modell der realen Produktion. Eine Anlage weist 
normalerweise keinen Eigenbedarf auf, das heiSt sie hat 
keinen Zielbereich, der erreicht werden mu£. Das Angebot 
der Anlage Si(t) ergibt sich zu: 

30 

Sift; - pi(t) 

Nun kann das Transport problem, wie es in dem vorausgegan- 
genen Abschnitt gezeigt ist, definiert werden. In FIG. 5 
35 mu£ lediglich pi ft; durch Si(t) und o t (t) durch d t (t) er- 
setzt werden. 

13 
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VON DEM VERTEILUNGS PROBLEM ZU DEM TRANS PORT PROBLEM IN DEM 
ALLGEMEINEN FALL. 

Die Modelle zum Projizieren des bisher gezeigten Trans- 
portproblems sind auf wesentliche Werte begrenzt. Es wird 
nun ein Modell erstellt, bei dem Sicherheitswarenbestande 
und Auftrage mit mehreren Prioritatsklassen berucksichtigt 
werden : 

t+R(t) 

Li(t) = m j rt+i?ct);+ E lof(T) 

r=t+l psPprognostiziert 
(Verallgemelnerung fur nichtganzzahlige R(t) siehe oben) 

Anmerkung: Wenn auch Durchgangsf rachten trjft; (festgeleg- 
te Anlieferungen) berucksichtigt werden, dann reduziert 
sich der Bedarf zu: 

di(t) = Oi(t) - tri(t) + Li(t) - Li(t-l). 

Immer dann, wenn die Durchgangsf racht gro£ genug ist, kann 
dieser Bedarf moglicherweise negativ werden, d.h. dieser 
Knoten verandert sich zu einem Angebotsknoten mit positi- 
vem Angebot : 

Si(t ) = -r Oi(t). - tri(t) + Li(t) - Li(t-l) 

Das bedeutet, da£ immer dann, wenn die festgelegte Anlie- 
ferung tri(t) grofi genug ist, ein Verteilzentrum zu einem 
bestimmten Zeitpunkt t ein Angebotsknoten (anstelle eines 
Bedarf sknotens) werden kann. In diesem Fall kann ein der- 
artiges Angebot eines Verteilzentrums i zu einem bestimm- 
ten Zeitpunkt t hinsichtlich der Transportzeit einem be- 
liebigen Bedarf eines anderen Knotens zugeordnet werden. 
Insbesondere die Transportzeit und die Transportkosten fur 
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die zukunftigen Bedarfe des gleichen Verteilzentrums i 
werden auf Null gesetzt. Zu Beginn des Planungshorizonts 
t = 0 kann der anfangliche Bestand von Verteilzentren 
durch festgelegte Anlief erungen tri(O) modelliert werden. 
FIG. 6 zeigt ein kleines Netzwerk rait einer Herstellungs- 
anlage 602 und zwei Verteilzentren 604, 606. In dem Ver- 
teilzentrum 1 beispielsweise liegt aufgrund eines sehr 
grofcen anfanglichen Bestands in Periode 1 ein Uberangebot 
vor. Dies kann nun dem Bedarf der Periode 2 608 oder dem 
Bedarf des Verteilzentrums 2 zugeordnet werden. 

UNTERSCHIEDLICHE TRANSPORTM&GLICHKEITEN 

In dem vorausgegangenen Abschnitt wurde ein sehr einf aches 
Netzwerk mit nur einem Transportmittel zwischen zwei Kno- 
ten berucksichtigt . Es wird nun die Veranderung unter- 
sucht, die zu dem Prozefc fiihrt, eine Losung zu finden, 
wenn verschiedene Transportverf ahren zur Verfiigung stehen. 
Verschiedene Transportverf ahren bedeutet nicht notwendi- 
gerweise verschiedene Transportmittel, es kann sich auch 
auf Alternativen bei dem Weg beziehen (beispielsweise fur 
den Fall von Lastwagentransporten) . Der langere Weg mufi es 
wert sein, das heifit er ist nur dann sinnvoll, wenn da- 
durch niedrigere Transportkosten entstehen. 

FIG. 7 veranschaulicht ein Distributionsnetzwerk mit 
mehreren moglichen Mitteln fur den Transport wie bei- 
spielsweise durch Lastwagen 702 oder durch Flugzeuge 704. 
Verschiedene Bogen entsprechen in diesem Fall verschiede- 
nen Transportmitteln . Die Strecken sind immer mit den ver- 
wendeten Transportmitteln und mit der Zeit 706 des Trans- 
ports (in Perioden) markiert. In der Darstellung des opti- 
mierenden Netzwerks korreliert jeder Bogen mit einer maxi- 
malen Transportdauer . Beispielsweise kann fur die Zuord- 
nung eines Angebots aus Periode 1 zu einem Bedarf in Peri- 
ode 2 lediglich ein Transportmittel gewahlt werden, das 
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zum Transport nicht mehr als eine Periode bendtigt. AuSer- 
dem ist es sinnvoll, nicht lediglich eine durchfiihrbare 
Moglichkeit, sondern die beste Moglichkeit zu wahlen. Man 
kann dann die grundlegende Annahme formulieren, dafi fur 
5 den Fall, dafc ein langsameres Transportmittel gewahlt 
wird, es auch das preiswertere ist. Diese zusatzliche An- 
nahme verandert im Prinzip deshalb nur die Kosten der Zu- 
ordnung in dem entsprechenden Transportproblem. Die Ver- 
bindung mit dem Transportproblem ist in FIG. 8 gezeigt, in 
10 welcher ein Optimierungsnetz dargestellt ist, in dem meh- 
rere Transportarten vorliegen. Die Projizierung der Trans- 
portkosten wird durch die verschiedenen genutzten Trans- 
portarten erweitert. 



15 TRANSPORT UBER ZWISCHENLAGEREINRICHTUNGEN . 

Im Weiteren wird zunachst der Fall mit nur einer Moglich- 
keit des Transports fur, jede direkte. Strecke zwischen zwei 
Knoten betrachtet. AnschlieSend wird die Veranderung be- 
riicksichtigt, die sich ergibt, wenn fur eine derartige 

20 Strecke mehrere Alternativen zugelassen werden. 

ZWISCHENLAGEREINRICHTUNGEN IM FALL EINER MOGLICHKEIT 

PES TRANSPORTS PRO STRECKE, 
Bisher wurde lediglich ein zweistufiges Distributionsnetz- 
25 werk berucksichtigt . Die Verbraucher (in diesem Fall die 
Verteilzentren) werden von den Herstellungsanlagen direkt 
belief ert. Oftmals benutzen Distributionsnetzwerke auch 
zentrale Zwischenlagereinrichtungen. Durch Einfuhren die- 
ser Einrichtungen konnen die Transportwege zu den Verteil- 
30 zentren nun teilweise zusammengebundelt werden. Anliefe- 
rungen konnen ebenfalls gebiindelt und somit moglicherweise 
preiswerter gemacht werden. Das zentrale Lagerhaus funk- 
tioniert wie eine Durchgangs station, urn die Anlief erungs- 
optimierung zu unterstiitzen. Aus diesem Grund kann die 
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Produktion der Anlagen den Verteilzentren direkt zugeord- 
net werden. Fur die Dauer des Gesamttransports (uber die 
. Zwischenlagereinrichtung) kann die Summe der Transportzeit 
zu dem zentralen Lagerhaus und von dort zu dem Verteilzen- 
5 trum verwendet werden. Die Transportzeiten schlieSen gege- 
benenfalls bereits die Ein- und Auslagerzeiten mit ein. Da 
das erf indungsgemafce Verfahren den Bereich der Abdeckungs- 
profile der Lagerhauser maximieren soil, sind moglichst 
fruhe Anlieferungen und deshalb kurze Transportzeiten von 

10 Vorteil. Aus diesem Grund wird immer der schnellste Weg 
zwischen zwei Knoten bevorzugt, selbst wenn er nicht not- . 
wendigerweise der direkte Weg ist. In der graphischen Dar- 
stellung der Transportzeiten ist f urn die kiirzeste Ge- 
samttransportdauer zwischen zwei beliebig gewahlten Knoten 

15 zu erhalten, die transitive Hulle der graphischen Darstel- 
lung zu berechnen. FIG. 9 zeigt ein Beispiel mit der kur- 
zesten Transportzeit zwischen zwei Knoten, oben ist das 
gegebene Netz dargestellt, unten die dazugehorende transi- 
tive Hulle. 

20 In dem ersten Flufcdiagramm 902 sind die verschiedenen 

Transportstrecken und -zeiten gezeigt. Es dauert drei Pe- 
rioden, um von Knoten 904 zu Knoten 906 zu kommen. Urn al- 
lerdings von Knoten 904 durch Knoten 908 zu Knoten 906 zu 
kommen, dauert es nur zwei Perioden. Unter Berucksichti- 

25 gung dieser Angaben wird das Flufcdiagramm 910 erstellt, 
das die Mindestzeiten zeigt, die erforderlich sind, um 
sich zwischen den Knoten zu bewegen. 

Wenn die Zeitdauer des schnellsten Wegs auf diese 
Weise berechnet worden ist, werden dann die damit verbun- 

30 denen Kosten berechnet. Analog zu den Zeiten setzen sie 
sich aus der Summe der Kosten fiir jeden Teil des Wegs zu- 
sammen. Gibt es mehrere schnelle Wege, werden die Kosten 
des ef f izientesten und/oder giinstigsten unter ihnen ver- 
wendet . 



17 




WO 99/67729 PCT/EP99/04229 

ZWISCHENSPEICHEREINRICHTUNGEN IM FALL VON ALTERNATIVES 
MOGLICHKEITEN PES TRANSPORTS FUR JEDE STRECKE . 
Wie oben erlautert, konnen mit dem erf indungsgemafcen Ver- 
fahren auch alternative B6gen in dem Distributionsnetzwerk 
5 beriicksichtigt werden. Beim Transport uber Zwischen3cnoten 
stellt sich die Frage, welcher Bogen fur die verschiedenen 
Segmente verwendet werden soil. Wenn alle Kombinationen 
zugelassen werden, dann ergeben sich bei mehreren Zwi- 
schenknoten schnell sehr viele Moglichkeiten des Trans- 

10 ports zwischen den Endknoten. Bei naherer Untersuchung al- 
lerdings stellt sich heraus, dafi nicht jede dieser Bogen- 
kombinationen sinnvoll ist. Zum Beispiel wiirde man norma- 
lerweise nicht in einem ersten Segment eine Exprefcverbin- 
dung und in einem zweiten Segment eine langsame Verbindung 

15 benutzen. Aus diesem Grund ist zur Eingabe der Verbindun- 
gen folgendes Verfahren gewahlt worden: 

Zunachst wird die Anzahl der maximal zulassigen Transport - 
arten in dem Netz spezif iziert . Die Transportarten ent- 
sprechen hier verschiedenen Geschwindigkeitsklassen (bei- 

20 spielsweise Normal und Exprefc) . Wenn nun zwischen zwei 
Knoten eine direkte Verbindung definiert werden soil, dann 
gibt man fiir jede Transportart (d.h. Geschwindigkeit , 
Klasse, usw. ) die Kosten und die Dauer an. Der Optimierer 
berechnet dann fiir jede Klasse getrennt den besten Weg, 

25 wie in dem letzten Abschnitt beschrieben. Diese Daten wer- 
den dann in dem Optimierungsnetzwerk verwendet. Bei Wegen 
mit mehreren Segmenten werden die Transportarten nicht 
miteinander vermischt. Wenn dennoch auf einer Strecke nur 
eine Transportmoglichkeit zugelassen werden soli, dann 

30 werden fiir alle Transportarten die gleichen Daten verwen- 
det. FIG. 10 zeigt ein Beispiel. Die Zeichnung zeigt oben 
das eingegebene Netzwerk, und unten das daraus durch den 
Optimierer erstellte Netzwerk. Der Weg 1002 zeigt ein Her- 
s t el lungs zentrum 1004, ein zentrales Lagerhaus 1006 und 
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ein Verteilzentrum 1008. Er zeigt einen Transportweg 1010 
per Zug, der 1 Periode benotigt, Transportweg 1012 per 
Lastwagen, der 3 Perioden benotigt, Transportwege 1014, 
1016 per Flugzeug, die jeweils 1 Periode benotigen. Die 
Bogen des Normalmodus sind mit Punkten gezeigt, die Bogen 
der Exprefcbetriebsarten sind mit Linien gezeigt. Zwischen 
den Knoten 1006 und 1008 gibt es nur eine reale Verbin- 
dung, namlich den Transport per Flugzeug. Der Weg 1018 
zeigt die von dem Optimierer zwischen dem Herstellungszen- 
trum 1004, dem zentralen Lagerhaus 1006, und dem Verteil- 
zentrum 1008 komprimierten Wege. 

UNZULASSIGE PERIODEN 

Im folgenden wird die Implement ierung von unzulassigen B6- 
15 gen erdrtert. Jede Transportverbindung und jeder Distribu- 
tionsknoten kann mit einer Liste unzulassiger Perioden 
versehen werden. Unzulassige Perioden sind Zeitfenster wie 
beispielsweise Feiertage, an denen der Transport runt oder 
an denen keine Waren angenommen oder ausgegeben werden 
20 konnen. Es wird an den oben beschriebenen klassischen Fall 
erinnert, in dem jeder Knoten in dem Optimierungsnetz fur 
einen Ausgangspunkt i, einen Startzeitpunkt t x , einen 
Zielort j und eine bestimmte Zielzeit tj einer bestimmten 
Materialzuweisung steht. Urn Feiertage zu beriicksichtigen, 
25 wird jeder Bogen soweit in die Zukunft geschoben, bis ti 
auf eine zul&ssige Periode zur Ausgabe von Waren beim Kno- 
ten i fallt (das heiEt, sowohl zu der Startzeit als auch 
zu der Zielzeit wird ein gemeinsames Zeitdelta addiert) . 
Danach wird tj ebenfalls solange verschoben, bis die Sum- 
30 men der zulassigen Transportperioden zwischen ti und tj fur 
den entsprechenden Transport ausreichen, und bis tj auf 
eine zulassige Periode zur Annahme von Waren beim Knoten j 
fallt. FIG. 11 zeigt ein optima erendes Netzwerk mit einem 
Horizont aus sechs Perioden 1102 , 1104, 1106, 1108, 1110, 

19 
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1112. Die Knoten sind mit der Distributionsnetzknotennuin- 
mer i und der zugehorigen Periodennummer ti markiert. Die 
unzulassigen Perioden sind jeweils schraf f iert/grau darge- 
stellt. In Periode 2 ist die Warenannahme im Distributi- 

5 onsnetzknoten 1 verboten. Die Warenabgabe ist in Knoten 2 
in Periode 5 unzulassig. In Periode 3 1114 ist der Trans- 
port unzulassig. Die beschriebene Distribution der unzu- 
lassigen Perioden ist naturlich ein Extremf all . In den 
meisten Fallen wird der Feiertag wenigstens die gleichzei- 

10 tig warenempfangenden und warenausgebenden Bereiche tan- 
gieren. 

DIE SUCHE NACH DEM MAXIMAL EN BEREICHSPROFIL . 

Das Ziel der vorliegenden Erfindung besteht darin, die 

15 Verteilzentren der Bedarf sknoten mit dem zur Verfugung 
stehenden Material derart zu belief ern, daE das Reichwei- 
tenprofil des Bestands uber alle Lager maximiert wird. Die 
Aufgabe besteht deshalb darin, eine geeignete Zuordnung 
von den Angebotsknoten zu den Bedarf sknoten zu finden. 

20 Bisher wurde dargestellt, wie eine Zuordnung zu einem ge- 
gebenen Bereichsprof il bestimmt wird, welche die oben an- 
gefuhrten Bedingungen erfiillt (minimale Transportkosten, 
die Feiertage usw. beriicksichtigen) . Dies geschieht durch 
die Abbildung eines Transportproblems und Losung desselben 

25 unter Verwendung eines effizienten Transportalgorithmus, 
beispielsweise des MODI-Algorithmus, der als Grundlage fiir 
den von der vorliegenden Erfindung vorgeschlagenen verbes- 
serten Algorithmus dienen kann. Nachfolgend wird das Pro- 
blem der Maximierung des Bereichs der Abdeckung definiert. 

30 
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MAXIMIERUNG PER REICHWEITENPROFILE BE I UBERANGEBOT 
Ein Uberangebot bedeutet, da£ dem geplanten Kundenbedarf 
vollstandig entsprochen werden kann. Somit hat der Trans- 
portoptimierer bei dem minimalen Reichweitenprof il von ei- 
5 ner Periode iiber alle Verteilzentren und Zeitpunkte hinweg 
eine Losung gefunden. Nun muS das Reichweitenprof il maxi- 
miert werden. Auf den ersten Blick besteht die einzige Op- 
tion darin, mit den gegebenen Moglichkeiten verschiedene 
Profile auszuprobieren und allmihlich auf das Optimum hin- 

10 zuarbeiten, z. B. durch ein Trial-and-Error-Verfahren. An- 
nahmen hinsichtlich der existierenden Angebots- und Be- 
darf ssituation, die die allgemeine Giiltigkeit des Prozes- 
ses nicht iibermaSig einschranken, konnen allerdings die 
Suche stark beschleunigen. Man mu£ allerdings das Risiko 

15 akzeptieren, dafi man nicht immer eine optimale Losung fin- 
den kann. Vor einer ausfiihrlicheren Erorterung sollte die 
Frage beantwortet werden, was unter einem maximalen Reich- 
weitenprof il zu verstehen ist. 

20 WIE SIEHT EIN MAX I MALES REICHWEITENPROF I L AUS? 

Wenn man die Reichweite aller Verteilzentren in dem Netz- 
werk beschreiben will, dann ist das Reichweitenprof il je- 
des einzelnen Verteilzentrums zu beschreiben. Da die er- 
f indungsgemafie Aufgabe allerdings darin besteht, die 

25 Reichweiten zwischen den L&gern zu einem gegebenen be- 
stimmten Zeitpunkt soweit wie moglich auszugleichen, 
reicht es aus, nur von einem globalen Reichweitenprof il zu 
sprechen, das dann fiir jedes Verteilzentrum als unterer 
Grenzwert giiltig ist. Genauer gesagt ist das globale Pro- 

30 fil R(t) das Minimum fiir das ortliche Profil Ri(t) uber 
den Knoten i. Die ortlichen Profile konnen moglicherweise 
zuweilen iiber dem globalen Profil liegen; sie konnen al- 
lerdings niemals darunter liegen, d. h. es gilt: R(t) = 
mini (Rt(t)}. 
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Wie bisher, wird im Weiteren nur von dem Reichweiten- 
profil gesprochen (das alle Verteilzentren betrifft) . Man 
wurde intuitiv vermuten, daS das maximale Bereichsprof il 
dasjenige Profil ist, das in jeder Periode maximal ist. 
Dieses trifft nicht zu. Die maximale Reichweite in einer 
bestimmten Periode karm sich signifikant von der Reichwei- 
te einer anderen vorausgehenden Periode unterscheiden. 

FIG. 12 zeigt als Beispiel ein Netzwerk mit zwei Her- 
stellungsanlagen und zwei Verteilzentren. Auf den Bogen 
sind die entsprechenden Transportzeiten eingetragen. In 
FIG. 13 sind die zugehorigen Produktions- und Bedarfspro- 
file gezeigt. Das Beispiel ist gem&£ dem einfachen Modell 
eines Distributionsnetzwerkes aufgebaut und beinhaltet 
keine Sicherheitswarenbestande , Initiallagerbes tande oder 
feste Zugange. Alle Bedarf sklassen weisen gleiche Priori - 
tat auf, und es gibt nur eine Transportart . Fur die fol- 
"gende Betrachtung des Reichweitenprof ils wird immer dieses 
einfache Modell verwendet. In einer bestimmten Periode ist 
die maximal mogliche globale Reichweite diejenige Reich- 
weite, die beide Verteilzentren in dieser Periode gleich- 
zeitig erreichen konnen. In Periode 2 weist die Reichweite 
einen Wert Drei auf. Dieser Wert wird nur dann erreicht, 
wenn die Produktion der ersten Periode von Werk 1202(40 
Einheiten) den Verteilzentren zu gleichen Teilen zugeord- 
net wird. Das Werk 1202 hat keinen EinfluE auf den Bereich 
in der zwei ten Periode, da das Material von diesem Werk 
1202 nicht vor der dritten Periode angeliefert werden 
kann. Bei dieser Zuordnung ergibt sich fur das Verteilzen- 
trum 1 in der vierten Periode automatisch eine Reichweite 
von Eins, da dieses Verteilzentrum von dem Werk 2 erst in 
der funften Periode erreicht wird und da die anfangliche 
Produktion der Herstellungsanlage 1 bereit "auf gebraucht ■ 
ist. FIG. 14 zeigt die maximale globale Reichweite fur 
diese erste Zuordnung. Die minimalen Lagerbes tande beider 
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Zentren 1206, 1208 zu Beginn jeder Periode sind in Klam- 
mern gezeigt. Es ergibt sich immer dann eine unendliche 
Reichweite, wenn der Lagerbestand. mindestens fiir alle fol- 
genden Perioden ausreicht, da fiir Perioden jenseits des 
5 Horizontes ein Bedarf von null Einheiten angenommen wird. 

Will man nun die globale Reichweite in Periode 4 von 

1 auf 2 anheben, bietet sich folgende Zuordnung an. Von 
den in Periode 1 von Werk 1202 produzierten 40 Einheiten 
werden 30 Einheiten dem Verteilzentrum 1206 und 10 Einhei- 

10 ten dem Verteilzentrum 1208 zugeordnet. Von den 40 Einhei- 
ten der ersten Periode in Anlage 1204 ordnet man minde- 
stens 20 Einheiten dem Verteilzentrum 1208 zu. Bei einer 
geeigneten Verteilung der restlichen Produktion von Anlage 
1202 (beispielsweise auf beide Zentren 1206, 1208' in glei- 

15 chen Teilen) wird die in FIG. 15 gezeigte globale Reich- 
weite erreicht . 

Diese Zuordnung hat den zusatzlichen Vorteil, dafc das 
Minimum des globalen Bereichs uber alle Perioden von 1 auf 

2 angestiegen ist. Dieses Beispiel zeigt, daS es nicht im- 
20 mer gunstig ist, bei Beginn die groStmogliche globale 

Reichweite aufzubauen. Dies verhindert manchmal die Mog- 
lichkeit, zu einem spateren Zeitpunkt ein Minimum der glo- 
balen Reichweite auszugleichen . Dieser ProzeE des Ausglei- 
chens des Bestands ist in jedem Fall wiinschenswert , urn ein 

25 gutes Gesamtprofil zu erhalten. Die Reichweitenprof ile in 
Fig. 14 und Fig. 15 miissen deshalb so eingestuft werden, 
dafc das Reichweitenprof il in Fig. 15 bei Bewertung sich 
als besser als das Profil in Fig. 14 herausstellt . Es 
reicht nicht aus, nur das globale Minimum der Profile zu 

30 vergleichen (wie man denken konnte) , da die anderen Peri- 
oden, bei denen eine Reichweite moglicherweise nur knapp 
uber dem Minimum liegt, bei einem Vergleich nicht beriick- 
sichtigt werden. Es wird fiir das globale Reichweitenprof il 
deshalb die folgende Ordnung definiert: 



23 




WO 99/67729 PCT/EP99/04229 



Einstufung fur das globale Bereichsprof il : 

Fur zwei globale Bereichsprof ile R 1 und R 2 gilt: 

5 R 1 > R 2 :<=> minfR 1 (t) \R l (t) *R 2 (t) } > minfR 2 (t)/R* (t) *R* (t) } 

Fiir die Vollstandigkeit dieser Einstufung ist zu zeigen, 
da£ wenn zwei Elemente nicht verglichen werden konnen 
(die Profile weisen bei verschiedenen Punkten das gleiche 
10 Minimum auf), es immer ein Element gibt, das grofier ist 
als die beiden anderen Elemente. 

DAS MONOTONE REICHWEITENPROFIL 

Es ist sehr zeitaufwendig, mit den bekannten Verfahren 

15 (Trial-and-Error-Verfahren) das optimale Reichweitenprof il 
zu bestimmen, wenn beispielsweise das folgende Verfahren 
angewandt wird: Beginnend mit einer minimalen gleichmafii- 
gen Reichweite von 1 zu jedem bestimmten Zeitpunkt vergro- 
fiert man die Reichweite in Periode 1 um ein kleines Delta 

20 der Reichweite. Wenn eine Losung gefunden wird, vergrofiert 
man die Reichweite in Periode 2 ebenfalls und so weiter. 
Nachdem alle Perioden behandelt worden sind, beginnt der 
Proze£ wieder bei der ersten Periode. Falls die Reichweite 
in einer bestimmten Periode nicht vergrofcert werden kann, 

25 wird diese Periode bei den weiteren Iterationen von den 

VergroSerungen ausgenommen . Die zeitliche Komplexitat die- 
ses Ansatzes wird durch 0(R max /R ain • H) Aufrufe des Trans- 
portoptimierers begrenzt, was moglicherweise zuviel Zeit 
in Anspruch nimmt, wobei R ma * als die optimale Reichweite 

30 definiert ist, der in der Periode auftreten kann; wobei 

Rmin als die Prazision bzw. Granularitat definiert ist, mit 
der dieser Bereich bestimmt werden soli, und H fiir den Ho- 
rizont steht. In diesem Fall mu£ der Optimierer immer die 
Daten des gesamten Horizonts optimieren. 
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Aus diesem Grund werden erf indungsgemafi zunachst 
monotone Reichweitenprof ile betrachtet. In einem 
5 zweiten Schritt wird dieser Ansatz auf nicht monotone 
Bereiche von Reichweitenprof ilen verbessert. Da bei 
einer monoton steigenden Funktion das Minimum immer am 
Anfang liegt, ist das Profil, das unter den monoton 
steigenden Prof ilen das Maximum ist, das erste (mit 

10 dem kleinsten t) , das einen grofieren Wert als die 
anderen Profile erreicht. Ein Transportplan, der auf 
der Grundlage des maximalen monotonen Bereichsprof ils 
berechnet wird, ist sehr gut, doch ist er leider nicht 
immer optimal. Dies zeigt ein Beispiel. Die Tabelle in 

15 FIG. 16 zeigt eine mogliche Situation in einem 
Distributionsnetzwerk mit einem Werk und zwei 
Verteilzentren. Die Transportzeit ist gleichmafcig eine 
Periode. 

Das maximale monotone Reichweitenprof il Ft* 1602 
20 weist in jeder Periode einen Wert von 1 auf (mit 
Ausnahme der letzten Periode, da fur den Bedarf von 
Perioden jenseits des Horizonts der Wert 0 angenommen 
wird) , da in Periode 5 kein grofcerer Bereich mdglich 
ist. Die Tabelle zeigt den Plan, der sich ergibt, wenn 
25 in dem Optimierungsnetz das Angebot dem Bedarf 
zugeordnet wird und wenn die Anlieferung sofort 
beginnt. Bei Z M 1604 beispielsweise wird die Produktion 
von 30 Einheiten in der ersten Periode den ersten drei 
Bedarf en des Verteilzentrums 1206 (s. Fig. 12) 
30 zugeordnet. Wenn alle 3 zugeordneten Anlief erungen 
sofort beginnen, d. h. in dieser ersten 
Produktionsperiode, dann findet ein Gesamttransport 
von 30 Einheiten zu dem Verteilzentrum 1206 statt. Der 
Transportplan 2? 1604 reicht fur das maximale monotone 
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Reichweitenprofil aus; er bildet aber eine 
tatsachliche globale Reichweite (das heiSt das Minimum 
der tatsachlichen Reichweite der Verteilzentren) von 
R' rt , 1606 , der kleiner ist als der maximal mogliche 

5 Bereich R 1608 (bei der Zuordnung Z) . Mit anderen 
Worten, der Plan Z M 1604 ist fur den maximalen 
monotonen Bereich sehr gut (er ist /nonoton-optimal) , 
ist aber nicht vollig optimal. 

Ein optimaler Transportplan fiir die erste Periode 

10 ergibt sich wie folgt: Die Zuordnung von Angebot der 
ersten Periode sollte keinen negativen Effekt auf den 
Bereich in samt lichen Perioden haben. Durch die 
folgende Annahme kann der Aufwand einer derartigen 
Suche gemaS diesem Plan wesentlich reduziert werden. 

15 Urn den monotonen optimalen Transportplan der ersten 
Periode zu erhalten, reicht es aus, die Suche nach dem 
maximalen Bereichsprof il auf solche Profile zu 
beschranken, die nur in dem Anf angsbereich nicht 
konstant sind, das heifct bei t < max ifj {Tij} mit Tij 

20 Transpprtzeit zwischen den Knoten i und j. Man mu£ das 
Reichweitenprofil nur bis zu der Periode der maximalen 
Transportzeit erhohen; in alien folgenden Perioden 
kann das Profil einen konstanten Wert aufweisen. 

Wenn bei der Bestimmung des 

25 Optimierungsnetzwerkes nicht der gesamte Horizont H 
betrachtet wird, sondern nur die Perioden t < bis 
zu einem bestimmten Zeitpunkt < H , das heifct 

wenn bei t <* lediglich Angebot st(t) und Bedarf 

di(t) beriicksichtigt werden, dann liefert das 
30 Optimierungsnetzwerk ein maximales Reichweitenprofil 
R 101 " , das in ieder Periode mindestens so groE ist wie 
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das maximale Profil liber den gesamten Horizont 

hinweg : 

Vt < : (t)> Rjoax (t> - 



SCHNELLE BEREICHSPROFILSUCHE 

Nachdem der kurze Horizont definiert worden ist, soli 
er auf eine Periode verkurzt werden. Das result ierende 
Optimierungsnetzwerk (das nur dann sinnvoll ist, wenn 
es moglich ist, den Transport in einer Periode 
auszufuhren) , ist sehr klein, und die Bogenf lusse 
konnen sehr schnell berechnet werden. Die maximale 
Reichweite in Periode 1 wird durch eine binare Suche 
ermittelt. In dem nachsten Schritt wird die zweite 
Periode hinzugefiigt, und es wird versucht, iiber beide 
Perioden hinweg ein maximales Bereichsprof il zu 
finden. Dabei wird das Reichweitenprof il aus dem 
ersten Schritt, das bis jetzt nur fur die erste 
Periode gilt, beibehalten, und die zweite Periode wird 
mit dem Wert der vorausgegangenen Periode (in diesem 
Fall der ersten) initialisiert Falls eine Losung fur 
dieses Profil gefunden wird, kann die Reichweite der 
zweiten Periode wieder unter Verwendung einer binaren 
Suche maximiert werden, usw. Falls fur einen kurzen 
Horizont die Reichweite R (H^ ) von R(H km -l) 

iibernommen wird, das heiSt falls keine Losung gefunden 
wird, wird die Reichweite in der Periode urn eine 

Periode gesenkt. Das resultierende Profil wird immer 
monoton gehalten, das heiiSt in alien Perioden t < 
wird die Reichweite auf R(t) = min{ R (t) , R ( } 

gesetzt. SpStestens jetzt muS eine Losung moglich 
sein, da sich eine Reichweite von 
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RiH^) = R(H btr7 -l) -1 aus dem Bereich RfH^ -1) in 
Periode H^z ergibt. Nun wird nochmals versucht, den 
Wert dieser und der folgenden Perioden zu erhohen, 
beginnend mit der ersten Periode, in der bei diesem 

.5 Schritt eine Reichweite verandert wurde. Der 
betrachtete Horizont bleibt bei . Nachdem = 

mzXfjfTyj} mit dem kurzen Horizont bestimmt wurde, 
werden im nachsten Schritt alle darauf folgenden 
Perioden wie eine einzige Periode behandelt. Der 

10 betrachtete Horizont wird auf den Gesamthorizont 
erweitert, und alle folgenden Perioden werden auf den 
Wert RiH^) gesetzt. Nach diesem Schritt ist das 
Verfahren beendet, und das maximale monotone Profil 
ist gefunden worden. Es steigt bis zu einem Punkt 7^ 

15 "an und bleibt danach konstant. 

FIG. 17 zeigt ein Beispiel gemaS dem 
erfindungsgemaSen Verfahren. Schritt 4 1702 zeigt 
einen Ruckschritt. Die Reichweite in Periode 
(Zeitintervall) 4 kann nicht auf dem Wert aus Periode 

20 3 gehalten werden und wird deshalb auf einen Wert von 
1,25 reduziert. Alle vorherigen Perioden mit einem 
grofieren Bereich werden auf diesen Bereich reduziert. 
Danach beginnt man wieder, den Bereich der betroffenen 
Periode zu vergroEern. In jedem Schritt stellt der in 

25 der letzten Periode des gerade betrachteten Horizonts 
fur die Reichweite gefundene Wert, das hei&t Rdi^), 
eine Obergrenze fur die Reichweite dar, die in dieser 
Periode erreicht werden kann. Da das maximale monotone 
Profil auch TransportplcLne toleriert, die es unmoglich 

30 machen, ein allgemeines maximales Profil 
(beispielsweise in diesem Abschnitt) zu erreichen, 
kann man diese Grenzwerte verwenden, urn einen 
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moglichen Fehler in dem berechneten monotonen Profil 
sichtbar zu machen und moglicherweise darauf 
interaktiv zu reagieren. 

1st Opt(H) die Laufzeit des Transportoptimierers fiir 
5 den Horizont H, dann ist die Laufzeit im 
ungiinstigsten Fall: 

OflogfR^/R^) - T m • Opt(H m )). 
R^ ist die maximale Reichweite, die in einer Periode 

auftreten kann, R^ ist die Genauigkeit 

10 ( Granular i t^t) , mit der diese Reichweite bestimmt 
werden kann, 7^ ist die maximale Transportzeit 
zwischen zwei Knoten (bei der langsamsten 
Transportart) und H total ist der Gesamthorizont (fiir die 
Gesamtlauf zeit einschlieSlich der Laufzeit fiir den 

15 Transportoptimierer) . 

Bei der praktischen Anwendung stellte es sich 
heraus, da£ ein Riicksetzschritt kaum erforderlich ist. 
Wenn man davon ausgeht, da£ Engpasse mit grofcerer 
Wahrscheinlichkeit am Anfang auftreten und dafi das 

20 Produktionsvolumen insgesamt das Bedarf svolumen 
geringfiigig iibersteigt (Push-Distribution) , dann 
findet in dem letzten Schritt (in dem alle Perioden 
nach r^in das Optimierungsnetz aufgenommen sind) kein 
Riicksetzschritt statt, und die Anzahl der Aufrufe des 

25 Transportoptimierers betragt nur log ( R m / ^ ) fiir den 
Gesamthorizont. Dann ergibt sich zusammen mit dem 
Aufwand fiir den Anf angsbereich folgende Laufzeit: 

OdogfR^/R^) • (T^ •Opt(T Bn )+Opt(H ulal ))) 

30 
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Das erf indungsgemafie Verfahren erzeugt schnell einen 
monotonen optimalen Transportplan fur die erste 
Periode. 

5 NICHT AUSRE ICHENDES ANGEBOT 

In dem vorausgegangenen Abschnitt wurde angenommen, 
daS die Produktion irraner ausreicht, um alle Bedarfe 
immer zeitig abzudecken. Dies ist natiirlich nicht 
immer der Fall . Wenn ein nicht ausreichendes Angebot 
10 vorliegt, muS die Bedarf ssituation auf geeignete Weise 
gelockert werden. Dies kann auf zweierlei Weise 
geschehen : 

1. der Bedarf selber wird reduziert, das heiSt einige 



wenige Kundenauf trage werden nur teilweise oder 



2 . Es wird eine Verzogerung bei der Erf ullung der 

Auf trage gestattet. 
Auch konnen beide Varianten kombiniert werden. Eine 
Kombination ist insbesondere dann sinnvoll, wenn ein 
20 nicht ausreichendes Angebot nicht lediglich ein 
kurzfristiges Problem ist, das heifit wenn die 
Produktion allgemein unter dem Bedarf liegt. In dieser 
Situation wurden sich Verzogerungen ohne 
Bedarf skiirzung weiterhin anhaufen. 



bedarfskOrzung 

Wenden wir uns nun dem ersten Fall zu (keine 
Verzogerung) . Um den Bedarf auf geeignete Weise zu 
30 verkurzen, werden Kundenauf tragen unterschiedliche 

Prioritaten zugewiesen. Eine Bestellung of (t) des 
Verteilnetzknotens i in der Periode t wird zusatzlich 
mit ihrer Prioritat P markiert. Die unwichtigen 
Prioritatsklassen (das heiSt Auftrage in diesen 
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uberhaupt nicht ausgefuhrt. 
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Klassen) werden soweit gestrichen, bis das zur 
Verfiigung stehende Materialvolumen ausreicht, die noch 
verbleibenden Auftrage rechtzeitig zu erfullen. Es 
wird der gleiche erf indungsgemafce Ansatz wie 
5 vorgenannt erwahnt. Der ganzzahlige Teil [R(t)J 
spezifiziert die Anzahl der Prioritatsklassen, deren 
Auftrage vollstandig erfullt werden. Der rationale 
Teil R(t)- l_R(t)J ist derjenige Teil, der von der 
nachsten Priori t&tsklasse erfullt wird. Zum Beispiel 
10 werden 4 Prioritatsklassen definiert: 

1. Phoch: Kundenbedarf e mit hoher Priori tat. 

2. Pnonnai : Kundenbedar f e mit normaler Priori tat. 
3- Pniedrg: Kundenbedarf e mit niedriger Priori tat. 

4. Pprognostiziert : Zusatzlicher Bedarf , der nur 
15 prognostiziert ist, aber noch nicht von einem Kunden 

bestellt worden ist 
Ein Abdeckungsbereich R(t) = 2,7 wiirde bedeuten, dafi 
zum Zeitpunkt t alle Bedarf e fiir Kunden mit hoher und 
normaler Priori tat erfiillbar sind. Ftir Kunden mit 
20 niedriger Prioritat betragt der Erfullungsgrad 70 %. 

Die erf indungsgemafie Definition eines maximalen 
Reichweitenprof ils garantiert, dafi kein Bedarf mit 
hoherer Prioritat zuungunsten eines Bedarfs mit 
niedrigerer Prioritat erfullt wird. 

25 

ANL I EFERUNG MIT VERZ OGERUNG 

Bei den obigen Uberlegungen war keine verzogerte 
Anlieferung gestattet. Dies bedeutet, da£ immer dann, 
30 wenn ein Auftrag nicht rechtzeitig angeliefert werden 
konnte, er uberhaupt nicht angeliefert wurde. 
Beispiel: Wenn ein Bedarf hoher Prioritat fiir t = 1 
und einen Bedarf niedriger Prioritat fur t = 2 
vorliegt und wenn das erste Angebot bei t = 2 ankommt, 
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dann wird der Bedarf mit niedriger Prioritat das ganze 
Angebot erhalten, aber der Bedarf mit hoher Prioritat 
.wird nichts erhalten. 

Dieses Problem wird erf indungsgemafi durch Einfuhren 
zusatzlicher Bogen mit reduzierter Transportzeit 
geldst. Die Anlieferung mit einer Verzogerungszeit A 
entspricht der Reduzierung der entsprechenden 
Transportzeit um A. 



Fur diese zusatzlichen Bogen werden die Kosten wie 
folgt definiert: 



15 



Strafe = A 2 * 2C max 



A : die Verzogerungszeit 

C max : maximalen Kosten der Bogen ohne 
Verzogerungszeit 
20 Der Ausdruck 2 C max garantiert, da£ durch einen 
Wechsel, d.h. durch Vertauschen der Zuordnung 

Si(ti)-> djftj) und Si-(ti-)-> dy{ty) 
durch Si(ti)-> dj'(tj-)und Si*(ti*)-> dj(tj), 

eine Verzogerung nicht verhindert werden kann. 

25 

Wenn es bei den vertauschten Zuordnungen keine 
Verzogerungen gibt, dann sind ihre Kosten < 2Cmax < A 2 + 
2Cmax und somit billiger als die Zuordnungen ohne eine 
Verzogerung. Der Optimierer wird aus diesem Grund 
30 niemals eine Losung mit einer Verzogerung produzieren, 
die durch einen Tausch verhindert werden konnte. 
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Analog dazu garantiert der Ausdruck A , da£ eine lange 
Verzogerung durch einen Tausch nicht verkiirzt werden 
kann. 

5 BESCHREIBUNG DER PROJEKTION PES DEPLOYMENTPROBLEMS IN 
DAS VERTE I LUNG S PROBLEM ■ 

Im weiteren wird die Darstellung der bisher 
hinsichtlich des Trail sportprobl ems gezeigten 
Verteilungsprobleme nochmals in mathematischer Form 
zusammengef a£t . 
FALL 1: UBERANGEBOT 
VORAUSGEHENDE BERECHNUNG 

Zuerst miissen Produktion und Bedarf des 
Optimierungsnetzwerks berechnet werden: 

Verteilzentren : 

Dj(t) = I o/rt; - trj(t) .+ Ljft-D-Ljft) } 
P 

Falls Djft) > 0 , dann gilt 
dt(t) = Dt(t) {Bedarfsknoten} 
Sonst gilt Sift) = -Djft) {Angebotsknoten} 

25 mit einem angestrebten Lagerbestand der 

Verteilzentren: 

LjftMl-R(t)-LR(t)J)W^ 1 

30 

L'jft,r)= mjft+r) + I Eof(r) 

r=t+l P 

Herstellnngsanlagen : 
35 s i (0)=p i (0)+L an/dnglich 

Vi>0:s i (t)=P i ft) 

33 
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DAS OPTIMIERUNGS PROBLEM 
Minimiere : 

5 .' 

" Transportkosten = 

> 

Unter den folgenden Nebenbedingungen : 

10 Beschranktes Angebot: V/ : s i (t i )=^y itjt 

»j ' 

Beschrankter Bedarf : V/ '; d t (t s )=2*y Hdtl 

u t > 

FALL 2A: VERRINGERUNG PES BEDARF S BE I UNTERANGEBOT 
15 VORAUSGEHENDE BERECHNUNG 

Zunachst mussen wieder Produktion und Bedarf des 
Optimierungsnetzes berechnet werden: 

Verteilzentren : 

20 

Dj(t) = Zoffl) - trj(t) + Lj(t-1)-Lj(t)} 

p&LR(t)J 

Falls Dfi) > 0 , dann gilt dtft) = Dfi) {Bedarfsknoten} 
25 Sonstgilt Si(t)=-Di(0 {Angebotsknoten} 



bei einem angestrebten Lagerbestand der 

Verteilzentren: 
30 Lj(tMl-R(t)-LR(t)J)Vfa^^^ L',(t/R(t) } 

t+r 

L)(t,r)= I Zof(t) 

t-f+/ P<PLR(t)J 

35 

Arunerkung: Der ganzzahlige Teil /r (t)J spezif iziert die 
Anzahl der Priori tStsklassen, deren AuftrSge 
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vollstandig erfiillt werden. Der rationale Teil R(t)~ 
LR(t)J ist derjenige Teil, der von der nachsten 
Prioritatsklasse erfiillt wird. 

5 Hers t el lungsanlagen : 

Si(0) = p^y+L"* 8 ** 
Vt>: Si (t) = P>(t) 

10 DAS OPT IMIERUNGS PROBLEM 
Minimiere : 
Transportkosten = 

15 Unter den f olgenden Nebenbedingungen : 
Beschranktes Angebot: Vi : 5 l Y/J = Xj>, l ^ 

Beschrankter Bedarf: Vy : d J (( i )-^y itJI 

,J i 

FALL 2B: VERZOGERUNG MIT UNTERANGEBOT 
20 VORAUSGEHENDE BERECHNUNG 

Wieder nuissen auf die gleiche Weise wie oben 
Produktion und Bedarf des Optimierungsnetzes zunachst 
berechnet werden, nur andert sich die 
25 Optimierungsfunktion durch Belegen der 

Verzogrungsbogen mit Strafkosten. Es konnen auch die 
zeitlichen Verzdgerungen durch Einfuhrung lediglich 
dieser zusatzlichen Bogen mit einer annehmbaren 
Verzdgerung begrenzt werden. 

30 
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DAS OPTIMIERUNGSPROBLEM 
Minimiere: 

. Kosten = Transportkosten + Verzogerungsstrafkosten 
Transportkosten = 



( 



\ 



5 



2 



Verzogerungsstrafkosten : 

£ £ (« * TiP«*-« "t j)' • 20^ . y M| 



10 Unter den folgenden Nebenbedingungen : 

Beschranktes Angebot: V/ ; s i (t i )^Yiy UJt 

JJ j ' ' 

Beschrankter Bedarf: V/ . d i (t i )-Yjy is 

J 

15 Maximalverzogerungsbeschrankung : 

Vi.tpj.t,: y^definiert <=> t s +T*~-- tjSA,, 

VARIABLE UND KONSTANTE 
Freie Entscheidungsvariablen : 
20 y ltJl = Teil des Angebot s von An 1 age oder Werk 

i zum Zeitpunkt der dem Bedarf des 
Verteilzentrums j zum Zeitpunkt tj 
zugeordnet ist 

25 Gebundene Variablen (mit beschrankter Entscheidung R(t) 

zu 

maximieren) : 

R(t) - Reichweite aller Verteilzentren zum 

Zeitpunkt / 

30 Lj(t) = Angestrebter Bestand im Verteilzentrum j 

36 
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sift) 

5 dj(t) 



zum Zeitpunkt t (als Funktion der 

Reichweite R(t)) 
= Angebot des Knotens (Produktionswerk) / 

zum Zeitpunkt / 
= Bedarf bei Knoten (Verteilzentrum) j 

zum Zeitpunkt t 



10 



Konstanten: 

tti(t) 



is l; 



anfdnglich 



= Produktionsrate in Produktionsanlage / 
zu einem bestimmten Zeitpunkt / 

= Zuvor festgelegtes Angebot (im 
Durchgang) im Verteilzentrum oder 
zentralen Lagerhaus j zum bestimmten 
Zeitpunkt / 

~-tri(0)= anfanglicher Bestand in 

Produktionsanlage i oder Verteilzentrum 



20 



o/ft) 



m/t) 



25 



= Bedarf spri or i tat P, Pe{Phoch>. • • ,Pprognosnzien}- 

{1, 2, 3, . . .} WObei Phoch~l< • * - < Pprognastiziert 

= Auftrag (=angehauf ter Bedarf) fur 
Prioritatsklasse P fur Verteilzentrum j 
zum Zeitpunkt t 

= Sicherheitslagerbestand im 
Verteilzentrum j zum Zeitpunkt / 



D 



= Transportart 
D e {langsamste, 



schnellste} 



30 Cjf 



= Transportkosten 
Materialeinheit 

37 



fiir Transport einer 
mit Transportart D 



WO 99/67729 



PCT/EP99/04229 



DschnelUte 
Dlangsamste 

langsamere(D) 



zwischen Knoten / und j 

Transport zeit fiir den Transport einer 
Materialeinheit mit Transportart D 
zwischen Knoten i und / 

Schnellste Transportart: VDJj: 

T DschnelUte < fp D 
I ij ij 

Langsamste Transportart: VDJj: T i j )Umgsamsu 
Nachstlangsamste Transportart nach 
Transportart D : 

-r3D ':[ VijiTf^^J 0 ^^] 

Maximal zulassige Verzogerung eines 

Bedarfs vom Verteilzentrum j zu einem 
bestimmten Zeitpunkt /, 



ALGORITHMEN 

Nachdem die Formulierung des Problems erfolgt ist, 
wird der anzuwendende Algorithmus erortert. Die 
grundlegende Idee des Algorithmus besteht darin, das 
Reichweitenprof il R(t) solange inkremental zu 
mutieren, bis das optimale Profil ermittelt ist. Die 
grundlegenden Verf ahrensschritte des Algorithmus 
umfassen die folgenden Punkte: 

1. Initialisieren das Reichweitenprof ils R(t) 
(z.B. R(t)=l). 

2. Berechnen der erf orderlichen Bedarfe dj (t) in 
den Bedarf sknoten, urn Prof ilbeschrankungen zu 
erf iillen. 

3 . Konstruieren des preiswertesten Flu£ fiir diese 
Bedarfe unter Verwendung eines beliebigen 
Minimalkostenf luSalgorithmus . 

4. Wenn keine Losung gefunden wird, Verkleinern des 
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Reichweitenprofils R(t), ansonsten VergroSern 
des Reichweitenprof il R(t) . 
5. Zuriickgehen zu Schritt 2, wenn die optimale 
Losung nicht gefunden wird. 

Schritt 2 ist der gleiche fur alle Algorithmen, wie in 

den Abschnitten oben beschrieben: 

dj(t) = Zof(t)-trj(t)+Lj(t-l)-Lj(t)} 
p 

wobei Lj(0=(l'R(t)-LR(Oj)'L)(tM^^ ) 
L'/t,r)= m/t+r) + £ lof(r) 

r=f+/ P 

Verschiedene erf indungsgemafce Alternativen zu Schritt 
4 werden in dem folgenden Abschnitt vorgeschlagen . 
Zunachst wird das Gesamtproblem erortert, d.h. die 
effiziente Konstruktion eines maximalen Bereichs von 
Reichweitenprof ilen. Dieser Algorithmus kann durch die 
Einschrankung auf einen monotonen Bereich von 
Reichweitenprof ilen, der unten beschrieben wird, 
betrachtlich beschleunigt werden. Im folgenden wird 
ein schneller Algorithmus fur monotone Profile fur die 
Maximierung des Reichweitenprofils angegeben. 

A. MAXIMALER BEREICH DES ABDECKUNGSPROFILS 

Zum Konstruieren des maximalen (ganzzahligen) Bereichs 
des Reichweitenprofils wird das tatsachliche Profil 
immer dann iterativ urn 1 vergrofcert, wenn dies fiir 
jeden Zeitschritt, der den gesamten Planungshorizont 
mehrmals uberstreicht , moglich ist. 

Algorithm us fiir Schritt 4 

While Prof ile_enlarged do begin 

{Den Planungshorizont uberstreichen} 
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enlarged = false; 



for all t (Planungshorizont) do begin 
R(t):=R(t)+ 1; i R vergrofcern } 

5 Generate min_cost_f low problem F(R) for profile R(t) 

{=2. Schritt (s.o. ) } 

If min_cost„flow (F(R)) solvable 
then enlarged = true 

else R(t):=R(t)-l {Rucksetzen}; 
10 . end; 
end; 

Die Verallgemeinerung zum Vergro&ern des Profils urn 
eine kleinere Schrittweite A < 1 anstatt 1 ergibt : 

15 

Alaorithmus von Schritt 4 mit Schrittweite A fu r 
Profil R 

While Prof ile_enlarged do begin 
20 {Den Planungshorizont iiberstreichen} 

enlarged = false; 



for all t < Planungshorizont do begin 
R(t):=R(t)+ A; {R vergrofiern} 

25 Generate min_cost_f low problem F(R) for profile R(t) 

{= 2 . Schritt (s.o. ) } 
if min_cost_flow (F(R)) solvable 
then enlarged = true 

else R(t):=R(t)-A {Rucksetzen} ; 

30 end; 
end; 
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w MAXTMALER MONOTONER BEREICH VON RKICHWEITEfl PROFILEN 

Die grundlegende Idee besteht darin, das Problem 
induktiv zu losen: 
5 • Initialisieren mit einem leeren Profil. 

• Konstruieren einer Losung fur t Zeitschritte 
und diese Losung fiir t + 2 Zeitschritte 
verbessern. 

• Bis t = Planungshorizont iterieren. 

10 

Alaorithrmis fiir S chritt 4 



R(t+l)=Rft); 

Generate min_cost_f low problem F(R) for profile R(t) 

{= 2. Schritt (s.o.)} 

If min_cost_flow (F(R)) solvable 

Then Begin {R(t + D 

vergrofiern} 



Initialize R (t) : 



Vt:R(t)=0 



15 



t-to; { Planungsstartzeitpunkt) 

While t < Planungshorizont do begin 

{Induktionsschritt} 



Repeat 



25 



R(t+l):=R(t+l)+l; 

Generate min_cost_f low problem F(R) for 
profile Rft) 



Until min_cost flow (F(R)) not solvable 




End 



30 



Else 



{R verkleinern} 



Repeat 

R^^Rft+l) 

t R = min{t | R(t) = /?™ - 1}\ 
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V?x R : R(t') = R maz -l; 

Generate min_cost_f low problem F(R) for 
profile R(t) 
Until min_cost_f low (F(R) ) solvable 

Aloorithmus fur Schritt 4 - monotoner Rgreich von 
Reichweitenorof ilen mi t Schrittweite A 

Die Verallgemeinerung zum Vergrofcern des Profils um einen 
kleineren Schritt A<1 anstatt 1 ergibt: 

Initialize R((): Vt:R(t)=0 

t=tO; {Planning_Start} 
While t < Planungshorizont do begin 
{ Induktionsschritt} 

Rfl+!)=R(t); 

Generate min_cost_f low problem F(R) for profile R(t) 

{- 2. Schritt 

(s.o.)} 

If min_cost_f low (F(R)) solvable 

Then Begin {Rft+D) vergroSern} 

Repeat 

R(t+l):=R(t+l)+A\ 

Generate min_cost_f low problem F(R) for profile 
R(t) 

Until min_cost_flow (F(R)) not solvable 

R(t+l):=R(t+l)-A {riicksetzen} 
End 

Else {R verkleinern} 

Repeat 

ST*: =R(t+l)\ 
t^minft'fRft'^R^-A} 
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W>t R : Rft'^R™ - A 

Generate min_cost_f low problem F(R) for 
profile R(t) 



Until min.cost^flow^fK;; solvable 



5 



C. FFFTZTENZVERBESSERUNG FUR ROLLENDES PLANUNGS SCHEMA 

Fur das rollende Planungsszenarium mufi lediglich bekannt 
sein, wohin die Produktion zum Startzeitpunkt t 0 , versandt 
werden soil, d.h. die Versandplanung y itQ j tj zum Zeitpunkt t 0 
10 sollte auf eine optimale Planung erweiterbar sein, 
allerdings ist die Versandplanung fur die folgenden 
Zeitschritte moglicherweise nicht optimal. Sie werden 
erf indungsgemafi durch Verschiebung des Planurigsf ensters 
korrigiert . 

15 im folgenden werden die Reichweitenprof ile auf den 

folgenden Typ eingeschrankt : GleichmaSig bis zur maximalen 
Transportzeit t max und konstant nach t max Solange dieses 
monotone Profil unlosbar ist, wird der letzte konstante 
Funktionsteil verringert. 

20 

Alaorithmus fu r fichritt 4 - rollendes Planunqg schema 
Vorgegeben ist die maximale Transportzeit, i.e. t raa x 



25 



Ermittle das maximale monotone Reichweitenprof il R max (t) 
innerhalb des Planungsf ensters tplanning_start , 
planning_start+ t ffla x] ; 



Initialisiere R(t) : {konstante Fortsetzung von R] 
Vt^M^R^ft) 



30 



W>t max Jl(t)=K m (t ma )\ 



B! m :^F! ma (t ma )\ 
verandern} 



{letzten konstanten Teil von R 
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Bestimme min_cost_f low problem F(R) fur das Prof il R(t) 
While min_cost_£lovi(F(R)) not solvable do begin 
t^minft/RCOZr*-!} 

R-=>R- - 1 

Bestimme min_cost_f low problem F(R) filr das Profil R(t) 

End; 

VF.RALLGEf ^-TWERUNGEN 

Zum Konstruieren von nicht ganzzahligen Kapazitatsprof ilen 
mufi die Schrittweite A der moglichen Werte fur den Bereich 
des Abdeckungsprofils R(t) begrenzt werden. Wenn die 
Schrittweite A klein ist, sind viele Iterationen 
erforderlich. bis das maximale Profil gefunden ist. In den 
folgenden Abschnitten wird beschrieben, wie die Suche nach 
maximalen monotonen Prof ilen durch einen binaren Suchprozefi 
beschleunigt werden kann. 

1 , pTTJ&RR STICHE WACH MAXI MRT.RM MONOTONEN PROFIftBfl 

Bei dem Algorithmus fur den maximalen monotonen Bereich von 
Abdeckungsprofilen mit Schrittweite A werden zwei Teile 
beriicksichtigt: 'R(t) vergroSern" and 'R(t) verkleinern" . 
In beiden Fallen wird die Suche nach der optimalen 
Abstimmung des letzten konstanten Teils des 
Kapazitatsprofils unter Verwendung einer binaren Suche 
beschleunigt. Eine einfache Erweiterung des 

erfindungsgemaEen Algorithmus mit einem binaren SuchprozeS 
wird nachfolgend angegeben. Aus Griinden der Deutlichkeit 
sind alternative Varianten ausgelassen (z.B. Wahl der 
Schrittweite S der binaren Suche mit durch Erfahrung 

gelernten Werten anstelle des kleinstmoglichen Werts 8= A) . 
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M^GORITHMPfi Ft TR SCHR TTT 4 - MONOTONER BEREJCH VON 
ARnECKUNGSPROFTLEN M T T SCHRITTWETTE A (BINAR5 SUCPE) 



While t < Planungshorizont do begin 
{ Indukt ionsschritt } 

10 Generate min_cost_f low problem F(R) for profile R(t) 

{=Schritt 2} 

If min_cost_flow (F(R)) solvable 
15 Then Begin {R(t+D vergroSern} 



5 Initialize R(t)\ 



Vt:R(t)=0 

t=to; { Planungsstartzeitpunkt } 



8:=A; 



Repeat 



{binare Suche-£ 



vergrolSern} 



R(t + 1) :=R(t+l)+6 



20 



5: =2*5) 

Generate min_cost_f low problem F(R) for profile 



R(t) 



Repeat 



Until min_cost_flow (F(R)) not solvable 

{binare Suche-£ 



25 verkleinern} 



R(t + 1) :=R(t + l)+8; 

5: =8/ 2; 



Generate min_cost_f low problem F(R) for profile 



R(t) 



30 



If min_cost_flow (F(R)) solvable 
Then 8: = l8i 
Else 8:=-l8l 
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Until /Sf=A; 

If min_cost_f low (F(R) ) not solvable 
' ** Then R(t+1) :=R(t + l)-8 

End 

5 Else ( R verkleinern} 

S:=A; 

Repeat {binare Suche-5 vergroSern} 

10 tf"*:=R(t+l) ; 

8:=2*8; 

Vt '2t R :R(f)-#" x -6; 
Generate min_cost_f low problem F (R) for profile 

15 R(t) 

Until min_cost_f low (F(R)) solvable 

Repeat {binare Suche-<5 

verkleinern} 
20 FP* x :=R(t + l) ; 

8: =5/2; 

Itxfrj&r ♦ 8/ 

Vt'ZtR-Rtt'i^F?*** 5; 

Generate min_cost_f low problem F(R) for profile 

25 J? ft; 

If min_cost_flow (F(R)) solvable 

Then S: = fsf 
Else 8:=-l8l 
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3. Wenn nicht das ganze Profil festgelegt wird, dann gehe 
zuriick zu Schritt 1 (und lose den restlichen Teil) . 

Fur Schritt 2 miissen die Zeitschritte t gewahlt werden, bei 
5 denen 



d.h. die stark mono ton ansteigenden Teile des maximalen 
monotonen Bereichs des Reichweitenprof ils konnen festgelegt 
werden. Aufcerdem kann auch das Ende des Planungshorizonts 
10 festgelegt werden: 



Deshalb kann das Planungsf enster mindestens urn 1 reduziert 
werden (eventuell um mehr, wenn wir auch gewisse 
festgelegte, stark monoton ansteigende Teile von 

15 reduzieren konnen) . Damit ist sichergestellt , da£ die obige 
Schleife hochstens T Schritte wiederholt wird, wobei 
T=planning_end-planning__start (Anzahl der Zeitschritte in 
dem Planungsf enster) , d.h. die zeitliche Komplexitat wird 
durch T*0(^ a) begrenzt, wobei 0(1?) die zeitliche Komplexitat 

20 des Algorithmus fur den maximalen monotonen Bereich des 
Abdeckungsprof ils ist. Die algorithmische Komplexitat 
steigt im Vergleich zu monotonen Profilen vorteilhaft 
hochstens um einen Faktor T , was es ermoglicht, auch 
umfangreiche Reichweitenprof ile mit herkommlichen 

25 Rechneranlagen iiberhaupt zu berechnen. 



F*(t) <FT(t+l) «> F?(t)=R(t), 



(planning_end) =R (planning^end) 
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Patentanspruche 



1. 

5 a. 
b. 
c . 

10 

15 



, Verfahren zum Maximieren der Lagerreichweite von 
Abdeckungsprof ilen, umfassend folgende Schritte: 
, Bereitstellen einer Lagerreichweitenfunktion fur das 

Abdeckungsprof il ; 
. Initialisieren der Lagerreichweitenfunktion mit 

mindestens einem Startwertparameter ; 
. Ermitteln einer optimalen Lagerbestandsmenge zur 
Losung der Lagerreichweitenfunktion fiir das 
Abdeckungsprof il, wobei das Berechnen folgende 
Schritte umfafit: 

i. Anwenden eines Minimalkostenf lu£algorithmus zum 
Bestimmen einer optimalen Transport losung zum 
Transport ieren der erf orderlichen Bestandsmenge ; 

ii. Bestimmen, ob eine Losung gefunden wurde, 
nachdem der Minimalkostenf lufcalgorithmus fiir den 
Bereich der Abdeckungsprof ilfunktion mit dem 
Startwert angewendet wurde; 

iii. inkrementales Verkleinern des Startwerts, wenn 
keine Losung bestimmt wurde, und Wiederholen des 
Berechnens einer Losung, bis eine optimale 
Losung fiir den Minimalkostenf lufialgorithmus 
ermittelt wurde; 

. inkrementales Vergrofiern des Startwerts, wenn eine 
Losung bestimmt wurde, und Wiederholen des 
Berechnens einer Losung, bis das Berechnen einer 
Losung zu einer optimalen Losung fiir den 
Minimalkostenf lu&algorithmus f iihrt . 

. Verfahren zum Maximieren eines Bereichs von 
Abdeckungsprof ilen nach Anspruch 1, weiterhin mit 
folgenden Schritten: 

. Ankoppeln an ein Online- 



48 



WO 99/67729 



PCT/EP99/04229 



Transaktionsverarbeitungssystem; 

b. Lesen von Transaktionsdaten von dem Online- 
Transaktionsverarbeitungs system; und 

c . Eingeben der Transaktionsdaten f ur den Bereich der 
Abdeckungsprof ilf unktion . 

3. Rechnerlesbares Medium mit mehreren darauf 
gespeicherten Anweisungen, wobei die mehreren 
Anweisungen Anweisungen en thai ten, die bei 
Ausfuhrung durch einen Prozessor bewirken, daS der 
Prozessor folgende Schritte durchfiihrt: 

a. Ankoppeln an ein Online- 
Transaktionsverarbeitungssystem; 

b. Lesen von Transaktionsdaten von dem Online- 
Transaktionsverarbeitungssystem; 

c. Eingeben der Transaktionsdaten in die Funktion fur 
einen Bereich des Abdeckungsprof i Is ; 

d. Initialisieren der Funktion fur einen Bereich des 
Abdeckungsprof ils mit einem Startwert; 

e. Berechnen einer erf orderlichen Bestandsmenge zum 
Erfiillen der Lagerreichweitenf unktion fur das 
Abdeckungsprof il ; 

f . Berechnen einer optimalen Ldsung, wobei das 
Berechnen folgende Schritte umf a£t : 

i. Anwenden eines Minimalkostenf luEalgorithmus zum 
Bestimmen einer optimalen Transport ldsung zum 
Transportieren der erf orderlichen Bestandsmenge; 

ii. Bestimmen, ob eine Losung gefunden wurde, 
nachdem der Minimalkostenf luEalgorithmus fur den 
Bereich der Abdeckungsprof ilfunktion mit dem 
Startwert angewendet wurde; 

iii. inkrementales Verkleinern des Startwerts, wenn 
keine Losung bestimmt wurde, und Wiederholen des 
Berechnens einer Losung, bis eine optimale 
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Losung ftir den Minimalkostenf lufcalgorithmus 
ermittelt wurde; 
g. inkrementales Vergrofiern des Startwerts, wenn eine 
Losung bestimmt wurde, und Wiederholen des 
5 Berechnens einer Losung, bis das Berechnen einer 
Losung zu einer optimalen Losung fur den 
Minimalkostenf luEalgorithmus f iihrt . 

4 . System zum Maximieren der Lagerreichweite von 
10 Abdeckungsprof ilen, umfassend: 

a. Mittel zum Ankoppeln an ein Online- 
Transaktionsverarbeitungs-system; 

b. Mittel zum Lesen von Transaktionsdaten von dem 
Onl ine-Transak tionsverarbei tungs sy s tern ; 

15 c. Mittel zum Eingeben der Transaktionsdaten in die 

Lagerreichweitenfunktion fur das Abdeckungsprof il ; 

d. Mittel zum Initialisieren der Lagerreichweite der 
Abdeckungsprof ilfunkt ion mit einem Startwert; 

e. Mittel zum Berechnen einer erf orderlichen 

20 Bestandsmenge zum Erfiillen der Funktion fur eine 

Lagerreichweite des Abdeckungsprof ils ; 

f . Mittel zum Bestimmen einer optimalen Losung, 
umfassend: 

i. Mittel zum Anwenden eines 

25 Minimalkostenf lufialgorithmus zum Bestimmen einer 

optimalen Transportldsung zum Transportieren der 
erf orderlichen Bestandsmenge; 

ii . Mittel zum Bestimmen, ob eine Losung gefunden 
wurde, nachdem der Minimalkostenf luSalgorithmus 

30 fur den Bereich der Abdeckungsprof ilfunktion mit 

dem Startwert angewendet wurde; 

iii. Mittel zum inkrementalen Verkleinern des 
Startwerts, wenn keine Losung bestimmt wurde, 
und Wiederholen des Berechnens einer Losung, bis 
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eine optimale Losung fiir den 

Minimalkostenf luSalgorithmus ermittelt wurde; 
g. Mittel zum inkrementalen VergroSern des Startwerts, 
wenn eine Losung bestimmt wurde, und Wiederholen des 
Berechnens einer Losung, bis das Berechnen einer 
Losung eine optimale Losung fiir den 
Minimalkostenf lufialgorithmus ergibt . 



51 



WO 99/67729 PCT/EP99/04229 

2 / 10 



FIG.3 




FIGA 




Lastwagen 



Flugzeug 



WO 99/67729 



6 / 10 



PCT/EP99/04229 



FIG.11 




FIG.12 

z-1206 



WO 99/67729 



7 / 10 



PCT/EP99/04229 



FIG.13 
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